A NUMERICAL APPROACH TO SOME BASIC THEOREMS IN 

SINGULARITY THEORY 



TA LE LOI AND PHAN PHIEN 

Abstract. In this paper, we give the explicit bounds for the data of objects 
involved in some basic theorems of Singularity theory: the Inverse, Implicit and 
Rank Theorems for Lipschitz mappings, Splitting Lemma and Morse Lemma, the 
bJQf density and openness of Morse functions. We expect that the results will make 

Singularities more applicable and will be useful for Numerical Analysis and some 
fields of computing. 



1. INTRODUCTION 

To make Singularity Theory more applicable it is important to make its basic 
results 'quantitative'. This direction of the theory is proposed by Y. Yomdin in 
[Yl] where he proves the quantitative Morse-Sard theorem by giving the notion of 
near-critical values of differentiable mappings and estimated these sets by the metric 
entropy. For the discussion on this direction and its developments we refer the readers 
to |Y2j and |Y-Cj and the references therein. We are interested in the numerical 
approach of this direction. In this paper, we give the quantitative versions of some 
basic theorems of Singularity theory: the Inverse, Implicit and Rank Theorems for 
Lipschitz mappings, Splitting Lemma and Morse Lemma, the density and openness 
of Morse functions. The explicit bounds for the data of the objects involved are 
estimated via the input data (e.g. C fe -norms, radii of the balls, ... ). The main tools 
that we use are some familiar methods of Singularities of differentiable mappings (see 
[A-G-Vj . |B-Lj . |G-G] , or [Ma] ), the quantitative forms of the Inverse and Implicit 
mappings theorems, and Morse-Sard theorem (see [Clj . |Paj . [P] , Theorems 13. II and 
13.51 in this paper, [Yl] and [Y-C] ). In our results, the estimates of the first order 
derivatives and the radii of the domains of the mappings involved are quite sharp. 
Since we use Lemmas 12.31 and I2.4[ the estimates of the higher order derivatives of 
the mappings involved are explicit but rather big. We expect that the results will 
make Singularities more applicable and will be useful for Numerical Analysis and 
some fields of computing. 

The plan of our paper is as follows: In Section 2 we recall some definitions and 
give the estimates of C fc -norms of compositions and inverses. In Section 3 we consider 
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the quantitative versions of the Inverse, Implicit and Rank Theorems for Lipschitz 
mappings. In Section 4 we give a quantitative form of diagonalization of matrix- 
valued mappings by upper triangular matrices, the quantitative versions of Splitting 
Lemma and Morse Lemma, and applications to the density and openness of Morse 
functions on a ball. 

2. PRELIMINARIES 

We give here some definitions, notations and results that will be used later. 
Let M mxn denote the vector space of real m x n matrices, 

||x|| = (l^i) 2 + • • • + |x n | 2 )^, where x G R", 

B"(xo) denotes the ball of radius r, centered at xq in R n , B™ = B™(0), and 
B" = B™, 

||A|| = max||a.|j =1 ||Ar||, where A G M mxn , or A is a linear mapping. 
Bmxn denotes the unit ball in M mxn , 

Sym(n) denotes the space of real symmetric n x n-matrices, 

A(n) denotes the vector space of all upper triangular n x n- matrices. 

Definition 2.1 (see |C2j ). Let / : W 1 — > M m be a Lipschitz mapping in a neighbor- 
hood U of Xq in M n , i.e. there exists a constant K > such that 



Then we denote L(f) < K. 

By Rademacher's theorem (see [F]), a Lipschitz mapping on a subset U of lR n is 
differentiable almost everywhere. The Jacobian matrix of the partial derivatives of 
/ at x, when it exists, is denoted by Jf(x). The generalized Jacobian of / at Xq, 
denoted by df(xo), is the convex hull of all matrices M of the form 



df(xo) is said to be of maximal rank if every M in df(xo) has the maximal rank. 
df(xo) is said to be of rank p if every M in df(xo) has rank p. 

Definition 2.2. Let /:{/—>• M. m be a differentiable mapping of class C fc , k > 1, on 
an open subset U of M n . Then the C fc -norm of / is defined by 



||/(x)-/(y)|| <K\\x-y\\, for all x,yeU. 



M=\im Jf{ Xi ), 



where (xi) converges to xo and / is differentiable at Xi for each i. 
For p < min(m, n), we denote 



dpxpf(xo) = {Mi G M pxp : there exists M 




\\f\\ ck = m a xsup\\DVf(x)\\. 

l<p<k -cf; 
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In the next sections we have to estimate the C fc -norm of compositions and inverses, 
to this aim we prepare the following two lemmas. 

Lemma 2.3. Let f : U — >■ V and g : V — > MP be differentiable mapping of class C k , 
k > 1, on open subsets U C M n , V C W m . Then 

\\9 o /||c* < (l k + 2 k + ■ ■ ■ + k k )\\g\\ ck maxdl/Uc*, ||/||^). 

In this paper, we denote 

E{Kf,K g ,k) = (l fc + 2 k + • • • + k k )K g m&x(K f , K k } ). 

Proof. By the Higher Order Chain Rule (see |A-M-Rj ). for p < k,we get the following 
estimation 

ii^o/)(x)n E -t^pw(^))iip ji /^)ii---ii^/wii 

st £ ^-rlMMI/lfe. 

<J> lbll^max(||/|| c ,,||/||^). 
i=i 

So ||^o /|| oh < (l fc + 2 k + ■ ■ ■ + k k )\\g\\ ck max(||/|| ofc , ° 

Lemma 2.4. Let tp : £7 — >■ V be a C k diffeomorphism between open subsets U, V of 
MJ 1 , k > 1. JTiere we /iai>e £/ie estimation 

y^Wck^EiiMc^WD^ik), 

where EI is constructed by the following recurrent method: 

Let M = E(\\ip\\ C k, max p\\\Dip- l \\ p+l , k - 1), Mi = \\Dip-\ 

0<p<k— 1 

and M p = S(Mp_i,M ,p- 1), forp = 2,--- ,k. Then EI(K, L, k) = M k . 
Proof. Let 

Inv: Gl(n) Gl(n) 

M ^ Inv(M) = M~ 1 . 

Since Dip^ 1 = Inv o £)<£> o <^ _1 , using Lemma |2.3[ we have the recurrent inequalities 
IIv^IIcp = niax(||L>v9 _1 ||, || Z^<^ — 1 1| op- 1 ) < ^(H^llc^- 1 , \\InvoD(p\\cp-i,p— 1), for p > 2. 

First, we estimate ||Invo£ty>|| *_i. From D p mv(M)(<5M) = p\(-l) p {M~ l 5M) p M- 1 , 
we get ||D p Inv(M)|| < p!||M _1 || p+1 . Therefore, using the notations K = \\(p\\ C k,L = 
||D<y9 _1 ||, we have 

\\D p Inv(Dtp(x))\\ < pllKD^x))" 1 !! 1 ^ 1 <p\L p+1 . 
Using Lemma [231 we get 

\\Inv o Dip\\ C k-i < E(K, max p\L p+ \ k - 1) = M . 

i<p<fc— i 
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Let 

Mi = L,M p = E(M p _ 1 ,M ,p-l),p = 2,--- , k. 
From || ip~ l \\ c i = H-Dc^ 1 !! < Mi, using Lemma [2.31 and recurrence, we have 

ll^llc* <M k = EI(K,L,k). 

□ 

Definition 2.5 (see |G-Lj ). Let L : M n — y M m be a linear mapping. Then there 
exist <Ji(L) > . . . > <J r (L) > 0, where r = rankL, so that L(B n ) is an r- dimensional 
ellipsoid of semi-axes (Ji(L) > ... > a r (L). Set a (L) = 1 and a r+ i(L) = ... = 
cr m (L) = 0, when r < m. We call (Jq(L), . . . , a m (L) the singular values of L. 

Remark 2.6. 

(i) ai(L) = \\L\\, a m (L) = min|| x [| = i \\Lx\\. 

(ii) If A G M is an eigenvalue of L, then a m (L) < |A| < o"i(L). 

Definition 2.7. A C fc function K on an open subset C/ of M™, k > 2, is 

called Morse if for every critical point x of /, i.e. Df(x) = 0, the Hessian Hf(x) 
is nondegenerate, i.e. a n (Hf(x)) > 0. 

3. THE INVERSE, IMPLICIT AND RANK THEOREMS FOR LIPSCHITZ MAPPINGS 

In this section, we present quantitative forms of the Inverse, Implicit and Rank 
Theorems for Lipschitz mappings, and give some explicit bounds in the smooth case. 

Theorem 3.1 (Inverse Mapping Theorem). Let f : M. n — y M n be a Lipschitz map- 
ping. Suppose that df(x ) is of maximal rank. Set 

S=l inf 1 



2 M iaf(x ) || Mo" 1 1| ' 
and r be chosen so that L(f) < K and 

df(x) C df(x ) + 5B nxn , when x G B"(x ). 
Then f : B^xo) — > f(B r p \( x o)) is a homeomorphism, where Pi — 2ic, fi^p^o)) 
contains B™ (f(x )), where p 2 = y, and L(f~ l ) < i on B™ (/(x )). 

Proof. See [Cl] and [F]. □ 

Remark 3.2. We make some comments on the pair (5, r) of the theorem that are 
often used latter. 

Let S = {A G M nxn : det A = 0}. Then by the Eckart- Young equality (see |G-L] ). 
we have 

— — = d(M, E), for every M in M nxn \ S. 
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Hence, 

5 = - inf - 1 1n = -d(df(x ), E). 
2 MoGdf(xo) \\M~ l \\ 2 v JK u; ' ; 

In words, 5 is half the distance from the generalized Jacobian of / at xq to the 

singular locus E. Note that if 5' < 5, then the theorem is also true when 5 is 

replaced by 5'. 

By the upper semicontinuous property of the generalized Jacobian (see [C2] ). for 
every 5 > there exists r > 0, such that 

df(x) C df(x ) + SB nxn , when x G B"(x ). 

So the quantity r reflects the rate of variation of the generalized Jacobian of / in 
a neighborhood of xq. If r' < r, then the theorem is also true when r is replaced by r' . 



Using Lemma [2.41 we have the following corollary. 

Corollary 3.3. With the assumptions and notations of Theorem l3.ll and in addition 
/ is a C k mapping, k > 2, and ^ K. Then we can choose 

Moreover f" 1 is also in class C k , < 4, and 

wr x \w<Ei{K,- 5 ,k). 

Remark 3.4. If F : U x V — > R™ be a Lipschitz mapping in a neighborhood of 
(xq, yo) in M. m x M n , then the generalized Jacobian of F at (xq, yo) satisfies 

dF{x , y )c{(M 1 M 2 ):M l e d x F{x , y ), M 2 e d 2 F{x , y )} , 

where diF(xo, yo) and d 2 F(xQ, yo) are the generalized Jacobians of F(-, yo) '■ U — > W 1 
and F(xq, •) : V — > M n at (x ,yo), respectively. 

Theorem 3.5 (Implicit Function Theorem). Let F : M. m x IR n — > ]R n 6e a Lipschitz 
mapping in a neighborhood of (x ,yo). Suppose that d 2 F(x ,yo) is of maximal rank 
and F(xq, yo) = 0. Set 

5 = 1 inf 



i : 
13 



2 Af 2 e9 2 F(xo,yo) (1 + (1 + iT) 2 || MP 1 1| 2 ): 

and r be chosen so that L{F) < K and 

dF(x,y) C dF(x ,y ) + 5B nx{m+n) , when (x,y) G B™ +n ((x , y )). 

Then there exists a Lipschitz mapping g : B p (xo) — > M n , where p = 2 (k+i) > an< ^ 
L(g) < y, such that 

g(xo) = yo, and F(x,g(x)) = 0, when x G B p (xq). 
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Proof, (c.f. jPaj.jP]). 

Set f(x, y) = (x, F(x, y)). Then / is Lipschitz with L(f) < K + 1 and 

df{x ,yo) C ( ( T J P.): Mi G diF(x Q ,y Q ), M 2 G a 2 F(x ,?/o) [> • 





M 2 ~ 



_ Mi M 2 t 

Since <9 2 .F(xo, yo) is of maximal rank, df(xo, yo) is of maximal rank. 

For M = ( Mi M 2 ) e d ^ x ^ we have M_1 = ( -M^Mi 
Therefore, 



i i 



IM" 1 !! = sup (\\x\\ 2 + \\M^ l M lX - M^ l y\\ 2 Y <{1 + {K + lf\\M^ l \\ 2 Y , 

||x-||2 + ||. y |P=l 



and hence 



Since 



1 1 
> 



M" 



(i + ^ + i^IIm^ 1 !! 2 ) 



i ■ 

2 



A = - inf n ; I \ n > 

2 Medf{x ,yo) IIM- 1 !! ~ 

by the supposition, we have 

df(x,y) C df(x ,y ) + A£ nxn , when (x, y) G B™ +n ((x , y )). 

By Theorem 13.11 / is locally invertible, with h(x,z)), (x,z) G 

B™ +m (x ,0), where p = ^ — ^ . Let y(x) = /i(x,0),x G B™(x ). Then g is Lips- 
chitz, g(x ) = y and F(x,g(x)) = 0. 
Moreover, when F is differentiable at (x,g(x)), we have 

\\Dg(x)\\ = ||-(f) _1 f(^^))ll 

< sup llMTl^ 

< sup ||M 2 _1 ||.K"< -if. 

n X n 

So L(y) < f . □ 

Corollary 3.6. With the assumptions and notations of Theorem 13 .51 and in addition 
F is a C k mapping, k > 2, and ||-F||c* fc < Then we can choose 

r,a„d r =A 



2(l + (l + ATIIf (0,0)-i|| 2 )* K 
Moreover g is also in class C k , \\Dg\\ < -j, and 

\\g\\ ck < C(K, k) = 2 k - 1 EI(K, f , k - l)K. 
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Proof. From F(x,g(x)) = 0, we have Dg(x) = - (jj^j ^(x,g(x)). 
By Theorem 13.51 \\g\lc 1 — \\Dg\\ < y. 

Applying the Higher Order Leibnitz Rule (see |A-M-Rj ) . we have 



\g\\ C k = max(||D^||, \\Dg\\ C k-i) 

%) llc*-i|l^-||c*-i 



^Ef *7 1 )ii(*r , iic"iiKi 



i=0 x ' 



(^) 1 lie- < ^(llf lie-, P (f )" - 1) < EI(K, f , fc - 1). 



To estimate || II c k - 1 ^ we use Lemma [2741 to get 

From this estimation, we get 

\\g\\ C H <2 k ~ 1 EI(K,f,k-l)K. 

□ 

Theorem 3.7 (Rank Theorem). Let f : M. n — > R m be a Lipschitz mapping in a 
neighborhood U of x G lR n with L(f) < K. Suppose that df(x) is of rank p for all 
x G U and that d pxp f(xo) is of rank p. Set 

5 = - inf 



2 Mied pxp f(x ) (1 + (1 + K) 2 1| M^ 1 1| 2 ) 2 
and r 6e chosen so that B™(xo) C J7 and 

<9/(x) C 9/(x ) + 5B nxn , for all x G B"(x ). 
Then there exist homeomorphisms 

<p : B»(x ) v(B^(rr )), and ^ : B£(/(x )) -> ^(B™(/(x )), 
w/iere pi = 2(j ^_ 1} , p 2 = y , < y and L(^) < 1 + f, 

such that ip(B n pi (x )) D B^ixo)), c^" 1 : B^(^(z )) -> R m , L(^ 1 ) < and 
if> o f o (p~\zx, . . . , z n ) = . . . , z p , 0, . . . , 0). 

Proof. Without loss of generality we can assume x = 0, /(0) = 0. 
Let tp : {R n , 0) -> (M n , 0) , be defined by 

<£>(2l, . . . , x n ) = (fi(x), . . . , fp(x), X p+ i, . . . , x n ). 

Then L(ip) < K + 1, and each M G cty?(0) has the form 

M =( M q M j 2 \ wh ere M x G 5 pXp /(0). 
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Hence d<p(0) is of maximal rank. By an estimation in the proof of Theorem 13.51 we 
have 



IM" 1 !! < (1 + {1 + K) 2 \\M{ 



1 1 1 2 > 



1 



1 1 1 

So S\ = - inf rr— , .. > S — - inf 

2 Med v (o) KM" 1 !! 2 Ahed pxp f(x ) (1 + (1 + K^HMf 1 !! 2 ) a 

inequality and the assumption, we have 

d<p(x) C dip(0) + 8 1 B nxn , for all x G B?. 



— r- From this 



Hence we can apply Theorem 13.11 to conclude that tp : B pi — > (p(B pi ) is a homeo- 
mor 
Let 



morphism, where p\ = 2 {k+i) > ^(BpJ D B P2 , where p 2 = y, and L((^ x ) < | 



2 ^ X )< 5 - 



^ = /o (/3 - 1 :B p2 ^M m . 

Then L(c/) < L(/)L(^ 1 ) < f , and g(z u . . . , z n ) = (z u ...,z p , g p+ i(z), . . . , g m (z)). 
Therefore, each M' G dg{z) is of the form 

M ' = ( k M 2 
Since (^(z) is of rank p when 2: G B p2 , 
(3.1) M 2 = 0. 

Let -0 : B p2 — > IR m , be given by 

/ 2/i \ / 2/i 



y P +i 



2/p+i - 9 P +i(yi, • • • ,2/ P ,0, . . . ,0) 



\ Vm / \ y m - g m (yi, ■ ■ ■ ,y P ,0, ■ ■ ■ ,0) / 

Then L(^) < 1 + < 1 + f , and each M" G <9V>(?/) is of the form 



M" 



I 

? I 



By Theorem 13.11 ^ is locally invertible. It is easy to see that ip is injective. So ip is 
a homeomorphism from B p2 onto its image. 

Because of (13. ip . gj(zi, • • • , -2 m ) — <7j(zi, . . . , z p , 0, . . . , 0) = 0, for every z G B P2 and 
j > p . Therefore ip o / o ip~ l = ip o g is represented by 

(z 1 ,...,z n ) h-> (z u ...,Zp,0, ...,0). 



□ 
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Corollary 3.8. With the assumptions and notations of Theorem 13 .71 and in addition 
/ is of class C k , k > 2, and ||/||c fc < K. Then we can choose 

andr K 



2(1 + (1 + ATIIWK)II 2 )* * 

K 



Moreover ip and ip are of class C , ||<£>||c* ^ K + 1, \\Dip\\ < -j + 1, and 

ck < C(K, 6, k) = E{EI{K + 1, -, k),K, k). 

o 



Proof. It suffices to estimate ||y?||c"= an d 1 1 V 7 1 1 o fc ■ 
By the definition of tp, we have ||^||c"= — 11/ ||c fe + 1 ^ K + 1. 
From < using Lemma 12. 4[ we have H^Hc* < EI(K +l,±,k). 

From this inequality, using Lemma \2.3\ we get 

IMIc* = 11/ ° y-lc* < + 1, k),K, k). 

Finally, by the definition of \\ip\\c k < IMIc* + 1> an d hence 

\\DiP\\ < \\Dg\\ + 1 < f + 1, and < E(EI(K + l,±,k),K, k) + 1. 

□ 

4. SPLITTING LEMMA - MORSE FUNCTIONS. 

In this section we give the quantitative versions of Splitting Lemma, Morse 
Lemma, the density and openness of Morse functions on a ball. 

To prove the Splitting Lemma, we prepare the following lemma, which gives a 
quantitative form of diagonalization of matrix-valued mappings by upper triangular 
matrices. 

Lemma 4.1. Let B : U — )■ Sym(n), be a C k matrix-valued mappings, k > 1, on a 
open neighborhood U of in W 1 . Suppose that \\B\\ C k — K, and 

B(0) = D = diag(±l,...,±l). 

Let 5 = — — — y= — — ; — . Then there exists a C k mapping 

4{K + l){K+2){l + {K + 2) 2 )n{n + l)) yy y 

Q:U S ^ A(n), where U s = U n B£, 
such that Q(0) = /„, B(x) = t Q(x)D Q(x), and 

\\DQ(x)\\ < (K+l){l + (K + 2)n(n + l)),\\Q\\ C k <C(K,n,k), 
where C(K, n, k) = 2 k ~\K + 1)EI(K +1,(K + 1)(1 + (K + 2)n(n + 1)), k - 1). 
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Proof. Let s = = dimA(n) = dim Sym(n). 

Consider F : U x A(n) ->■ Sym(n), F(x, Q) = B(x) - l QD Q Q. 
Then F eC k , F(0, I n ) = G M nxn , and 

BF BF 

— (0,/„) : A(n) Sym(n), —(Q,Q(H) = -*#D - D #- 
Denote # = (%), G = = -{ l HD Q + D„i?). Then 

Thus l§(0, I n ) is invertible, and 



HS(o,/»)ll = 

||ff|| = l ||tf||=l 



f|(0,/ n )||= sup tfHDo - D H\\ < sup 2||fl-||||D || = 2, 



ll(g(o,/»))- x || < ||(Sg(0,I„)H|* < = 

We are using Implicit Function Theorem, so we estimate some numbers related to 

F. From 

(4.1) 

F(x + Ax,Q + H) - F(x, Q) = B(x + Ax) - B(x) - 'HDqQ - 'QDqH - "HQDqH, 
we have 

BF — OF — — 

— (x, Q) = DB{x), —(x, Q)(H) = - t HB(x)Q - t QB(x)H. 
ox oQ 

For < r < v ^^ , when \\(x, Q)\\ < r, we have 

\\DF(x,Q)\\ < (\\DB(x)\\ 2 + sup \\ t HB(x)Q + t QB(x)H || 2 )i 

ll^ll=i 



< (K 2 + 2\\B(x)\\ 2 \\Q\\ 2 y <KVT+2r*<K + l. 

So L(F) < K + 1 on B" +s . To apply Theorem E31 to F, with 

1 



2(l + (l + (K + l) 2 s)) 



1 1 

2 



we chose r = min( ^~| X , = ^ to have 

||£F(x,Q)-£>F(0,I n )|| < (A-+l)r = 5i, when ||(x,Q)|| < r. 
According to Theorem 13.51 and Corollary 13. 6[ there exists a C k mapping 



Q:U s -> A(n) in class C\ U s = U n B 



5 ! 



where o = = — ■= r~r= — r= tt^, such that 

2(K + 2) 8(K + l)(K + 2)(l + (K + 2) 2 )s' 



Q(0) = I n , F(x,Q(x)) = B(x)- t Q(x)D Q(x) = 0, and \\DQ(x)\\ < '' ' 



St 
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Moreover, from (14. ip . we have 

Q)(Ax, AQ) = - 1 AQDqQ - t QD AQ, and hence ||g|| < 2||Q|| < 2, 
gfi(x,Q)(Ax,AQ) = -'AQDoAQ, and hence ||g(x,Q)|| < 1, 
(x, Q) = 0, when a G N» /3 G N s , a ^ ^ /? or \(3\ > 3. 

So ||F||(7fc < K + 2. Using Corollary 13.61 we get 

HQlIc* < C(K, n, k) = 2 k -\K + 2)EI{K + 2, fc - 1). 

□ 

Applying the above lemma and Implicit Function Theorem, we can get a quan- 
titative form of Splitting Lemma. 

Theorem 4.2 (Splitting Lemma). Let f : U — > R be a C k function on a neighbor- 
hood U of x in W 1 , k > 3. Suppose that ||/||c fe < K, and 

Df(x ) = 0, rankD 2 f(x ) = p. 

Let 

5/2 2CT 2 a 3 / 2 

5 = — ; — - — ---77- min(l, — — -, ; - p -), where a, D = crJD 2 f(x )). 

32(K+ I) 9 / 2 v '3(p 2 +p+l)'2(p 2 +p+l) ; ' p pK JK 0JJ 

Then there exists a C k ~ x diffeomorphism 

<p:B^ ^(BJ), with \\D<p\\ < 32( ^ 1)5 , 

Up' 

such that 

v 

f <p(x,y) = f(x ) + ^2±x 2 + a(y), x = (x u . . . , x p ) G W, y G 

i=i 

where a is of class C k and a(0), Da(0), D 2 a(0) vanish. 
Moreover, there exists a constant M(K, a p , k) > 0, such that 

| c *-i < M(K,a p ,k). 



bn-p 



Proof. We can assume x = 0, f(x ) = 0, and can choose the coordinate system 
(x,y) G 1R P x R q ,p + q = n, x = (xi, . . . ,x p ),y = (y x , . . . , y q ) so that 

/ d 2 f 

A = fri/(0,0)= —^-(0,0) 

is of rank p. Note that if cri >•••> cr p > are the singular values of A, then 
— >•••> — > are the singular values of A -1 and 

■i<iu-'ii<l. 



01 cr. 



v 
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Step 1. Consider the equation: ^(x,y) = 0. 

We have §£(0,0) = 0. To apply the Implicit Function Theorem 13.51 to §£ , we will 
determine some numbers. Let 

t= ! ^,a„dr< 8 



2(1 + (K + K 

P 

Then, for \\(x, y)\\ < r', 

\\D(^)(x,y)-D&(0,0)\\<Kr' = 8'. 
ox ox 

So we can apply Theorem 13.51 and its corollary 13.61 to have a C k ~ 1 mapping 



where 



r'5' 



such that 



2{K+1) SK(K+l)(l + (K + iy±Y 

p 

df 

5(0) = 0, ^(jg(y),y) = 0, 
and 

(4.2) \\Dg\\ < WgWc^ < 2 k - 2 EI(K, j-,k — 2)K = M X {K, a p ). 

Let a(y) = f(g(y),y)- Then a(0) = 0,Da(0) = 0, and 

df df df 

Da(y) = -K-(g{y), y)Dg(y) + —(g( y ), y ) = —(g( y ), y ). 
ox oy oy 

Since the mapping on the right side is of class C fc_1 , Da is of class C fc_1 , and hence 

a is of class C k . Moreover, %(g(y),y) = and 0(0, 0) = imply D 2 a{0) = 0. 

Let fi(x,y) = f(x,y) - a(y), (x,y) G B£. We have 

fMv),y) = 0, ^(g(y),y) = 0. 

Step 2. Let tpi(x,y) = (x + g(y),y), (x,y) G . Then ipi is a C k ~ x diffeomorphism 
from B^ to its image, and from (14. 2p . we get 

K 

Ti 

Let /2 = fx o if i. Then f 2 is of class C fc_1 , and 

/ a (0,y) = 0, ^(0,y) = 0. 

Note that 

y ) = ?Af x + ( ) ) an d ^ = i^i^A^i _ d 2 fi _ d 2 f 
dx dx ' &r 2 cte <9x 2 <9x cte 2 <9x 2 

Step 3. Let Qq G Gl(p) be the linear transformation so that 

t Q Q AQ Q = D Q = (±l,...,±l). 



(4.3) \\Dipx\\<l + \\Dg\\<l + —, \\<pi\\ck-i < 1 + Mi{K, a p ) = M 2 (K, a p ) 
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Moreover, choose Qo = SU, where U is an orthogonal matrix and S is a diagonal 
matrix, so that 

iiQoii 2 =-,iiQo 1 ir= °i=uw <k. 



Let B : B£ — > Sym(p), defined by 



1 ,1 Q2f, 



B(x,y) = {bij{x, y)) x< id < p , where % (x, y) = J j q x 'q x (stx,y)dsdt. 

Then B is of class C fc_1 , and 

f 2 (x, y) = 'xB(x, y)x, and B(0, 0) = A = HJ(0, 0). 

Set 

B(x,y) = t Q B(x,y)Q . 

Then B : B£ Sym(p) e C* -1 , 5(0,0) = D , and ||lXB|| < K = ^. According 
to Lemma [4.11 there exists a C fc_1 mapping 

G:Bj,->A(p), 

where 5 2 = min(<5i, — ■= — — — —= — — — ), 

such that 

Q(0) = I p ,B(x, y) = t Q(x)D Q(x), \\DQ\\ <(K + 1)(1 + (K + 2) 2 p(p + 1)), 

and 

(4.4) \\Q\\ c *-i <C(K,p,k-l) = M 3 {K,a p ). 

Let (p 2 (x,y) = (Q(x)QQ 1 x,y), (x,y)eB$ a . 

We are applying the Inverse Mapping Theorem 13.11 to Lp 2 ■ So we have to calculate 
to determine the pair (5, r) (see Remark 13 .2 j) and some numbers. First we have 

D<p 2 (x,y)(h,e) = (Q(x)Qo 1 h + DQ(x)hQ^ 1 x,e). 

Hence 

D<p 2 (0,0)(h,e) = (Q^Ke), and D<p 2 (0, 0)" 1 ^, e) = (Q H,e). 

Thus 



1 1 



\D<Pz(P, 0)- 1 !! sup (HQo/if + l)^ VIIQo|| 2 + l V 1 + ^ 

l|H||=l 
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So we get £ 3 = j- , / ° v . Let r 3 - 



2Vl + 0p 2(A + l)(l + (iT + 2) 2 p(p + l))VA 

Applying the Mean Value Theorem and (14.41) . when ||(x,y)|| < r3, we have 

\\Dtp 2 {x,y)- ^(0,0)11 = sup || Q(x)go^ + ^2(^)^0^ - 2(0)<3o ^11 

< IIQ(^)-Q(o))||||g 1 || + || J DQ(x)||||g 1 ||||x|| 

< 2||DQ||r 3 ||go 1 || 

< 2(K + l)(l + (K + 2) 2 p(P+±))VKr 3 = 5 3 . 

Now applying the Inverse Mapping Theorem 13.11 to </?2, we have 1 : B^ - — > K n . 
where 



(4.5) 5 = min( V 3)4 and ||^|| < I = 2 / + 



<5i y cr p 

Let 

V » = ¥ 7 1 o^ 1 :B?-> V (B?). 

Note that we used L(^ 2 ) = (A" + 1)(1 + (A + 2) 2 p(p + 1)) > 1, so, by Theorem EH 
y?2 1 (B^) C B£ (the domain of tpi). By the C k ~ 1 coordinate transformation (p, we 
have 

fo(p(x,y) = fifaiiv^fayfi+abj) 

= f 2 o(f2 1 (x,y) + a(y) 
v 

i=l 

Using (14. 3 p (14.51) and cr p < A, we can easily get the following estimate 



\\D<p\\ < H^ill \\D& 1 < (1 + 8A 2 (A + 1)(1 + ^))2 A /i±^ 

Op 

Moreover, using the Leibnitz Rule, we have 

Hv^llc*- 1 < IIQQo 1 !!^- 1 + 1 

<2 fc - 1 ||Q|| cfc - 1 ||g - 1 || + 1 

< 2 k - 1 M 3 {K,a p )VK + l = M 4 (K,a p ). 
From this estimation, (14. 3 j) (14. 5 p and using Lemmas 12. 3[ 12. 4[ we get 

ll^llc*- 1 = II ^2 1 ! I c fe — 1 

< EiW^Wo^-h ||^i||o*-i, *; - 1) 

< E(EI(M A , 2y^, k — 1), M 2 , k — 1) — M(K, a p , k). 

Step 4. To avoid the complicated formula for 5 we make some elementary estimates. 
Keeping track of the numbers during the proof, from 14.51 we have 

5 = ^5 3 min(a,b, c), 

where 
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1 



a 



2K(K + 1)(1 + 

a p 

1 



(^ + 1)(^ + 2)(1 + (^ + 2)Mp+1)' 

c ~ yr+^(^ + i)(i + (^ + 2)Mp + i))v / ^' 

Use cr p < K to get 

n - a \ > ^ 

" X(^ + 1)((72 + (X + 1) 2 ) ^ 4(if + l) 4 ' 

4 4 

u ' (K+a p )(K+2a p ){al+(K+2a p yip(p+l)) ^ 6(/ST+l) 4 (p 2 +p+l) ' 

473: 



c _ __^ = v J P U P > V"P"p 



So 



% /I+^( J fS'+ ( T p )(a2 + (/r+2 ( 7 p ) 2 p(p+l))v / X 2(E-+l) 4 (p 2 +p+l) ' 



2al a 3 J 2 



X > I /_3e ^£ minfl _£ l£ 

^ s V CTp +i 4(K+i) 4 mm ^ 3(p2 + p + !) ' 2 (p 2 + p+l 

rl 12 2al nl /2 



> — 7T7 77-7-7 mill ( 1 , 



p a v 



32(K + I) 9 / 2 v '3(p 2 +p+l)'2(p 2 + p+l) 7 ' 
We reduce the radius of the domain of ip to the last number to use in the statement 
of the theorem. □ 

Theorem 4.3 (Morse Lemma). With the assumptions and notations of Theorem 
1^-4 when p = n we have 

n 

f o ip( Xl , •••,&«) = f(x ) + ^2 ±x l- 

i=i 

Applying the quantitative Morse-Sard Theorem (see [Ylj or |Y-Cj ) and the In- 
verse Mapping Theorem, we give here a version for the density of Morse functions 
on a ball (c.f. |Y2l Th. 4.1, Th. 6.1]). 

Theorem 4.4. Fix k > 3. Let f : B n — > R be a C k -function with ||/o||c fc — K- 
Then for any given e > 0, we can find h with \\h\\ C k < e and the positive functions 
tpi, ^2, 1P3, d, M , N , n depending on K and e, such that f = f + h satisfies the 
following conditions: 

(i) At each critical point x\ off, o~ n (Hf(xi)) > ipi(K,e). 

(ii) For any two different critical points x% and Xj of f , \\xi — Xj\\ > d(K,e). 
Consequently, the number of critical points does not exceed N(K,e). 

(iii) For any two different critical points X{ and Xj of f , \f(xi) — f(xj)\ > if)2{K,e). 

(iv) For each critical point X{ of f , there exists a C k ~ 1 coordinate transformation 
ip : B$(xi) R n such that 

f ° • • • ,y n ) = vl + ■ ■ ■ + yf - yf+i vl + const i 



16 
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where 6 = ip^(K,e) and \\ip\\ck-i < M(K,e). 

(v) // ||D/(a;)|| < n(K,e), then x G B^Xj), with Xi is a critical point of f . 

Proof. In |L-P] , our proof of the theorem needs some corrections. Moreover, we can 
apply the Splitting Lemma [4.21 to get an alternative proof of (iv) in that paper with 
more explicit estimations for S and M. For these reasons, we make some improve- 
ments in detail in this present paper. 

(i) We are applying the results of Chapter 9 |Y-UJ to Df . For 7 > 0, denote 

7 = (Ai, A2, • • • , A n ) = (K, K, ■ ■ ■ ,7). Then, by definition, the set of 7-critical points 
and the set of 7-critical values of / are 

£(£>/„, 7, B*) = {x G IT : <Ti(D{Df )(x)) < \, i = 1, • • • , n} 

_ = {x G B" : a n {H f (x)) < 7}, and 
A(D/ ,7,fT) =/(E(D/ ,7,B n )). 

For a relatively compact subset A of ~R n , and r > 0, denoted by M(r, A) the minimal 

number of balls of radius r in IR n , covering A. 

Let e > 0. Applying Theorem 9.6 of |Y-Cj . when < r < s, 

M(r, A(D/ , 7, B") n B«) < c f ELo niin (a • • • A^ (^) 

n r / i_ \ / n 

[EE? mm (ffi (;)') + min (J) 

where A = 1, c = c(n, A;) and Rk(fo) — <k-iy. • If < r < 1 and r < Rk ^°i} e , then by 
taking the min and simplifying the right-hand side we get 

M(r, A(Df , 7, B") nB»)<c f £ ^(^( /o )|)n-_I_ + #n_L j . 
When 7 = r 1_ £, we have 

_ n 1 

(4.6) M^Ap/oj.FjnB^c^^W^r— 



i=0 ' 



Note that, by the definition of M(r, A), it is easy to see that M(2r, A r ) < M(r, A), 
where A r denotes the r-neighborhood of subset A of M. n . Therefore, if 

(4.7) M(r, A(Z)/ ,7,B n ) fl B™)m(B 2r ) < m(B"), 

where m(A) denotes the Lebesgue measure of A, then there exists vq G B™, such that 
t>o is not contained in a union of balls of radii < 2r that covers the r-neighborhood 
of A(Df , 7, tf) n B™, and hence B?(v ) n A(D/„, 7, B™) = 0. 
We want to find r, < r < min(e, 1, £°k £ ) satisfying 14.71 Combining 14.61 and 14. 7[ 
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we look for r satisfying 

" 1 F n 

Z-j v KWU7 ; r «-i+| (2r) n 

or 



^ V 2 ^Er=o^(^(/o)^) 

Taking 

r(K,e) =l m in(e,l,^^,( £ - T V^Y and 

7 (tf,e) =r(K,e) 1 -K 
we get 14.71 Then we can choose v G B™, such that B" (K £) (v) C B™ and every v' in 

2 

B " (Ke) (f ) is a 7 (if, e)-regular value of Df . 

Now, let / : R n — > R be a linear mapping with Dl = —v and fi = fo + I- Then 
< e- 1 ^, Afi = I>/ -v, and Hf x = Hf = D(Df ). So each v> € By,^) 

2 

is a 7(iT, e)-regular value of Dfi. In particular, at each critical point Xj of /i, we 
have 




(4.8) H^/iWII > <r n {Hf x (xi)) > -y(K,e). 

In other words, the smallest absolute value of the eigenvalues of the Hessian of fi 
at its critical points is at least ipi(K, e) = ^y(K,e). 

(ii) We are applying the Inverse Mapping Theorem 13.11 to Dfx : B™ — > W 1 at the 
critical points of f±. Let X{ be a critical point of f\. By (14. 8p we have 

' ' >\l{K,e). 



2||F/ 1 (x i )- 1 || - 2 
i? 



Choose 5' = hj(K, e), and r' = — . Applying the Mean value theorem, when \\x 



Xi\\ < r, we have 

WDiDf^-DiDf^W = \\D(Df )(x)-D(Df )( Xl )\\ <K\\x- Xi \\ < Kr' < 5' . 
Thus, by Theorem l3.ll Dj\ is invertible on B™, s , (xj) = B n 2 (xj). Hence, -D/f 1 (0) f~l 



7 2 (K,e) 

B" 2(Ke) (xj) has only one point, i.e. Xj is the unique critical point of /i in the ball 

SK 2 

B™ 2(if e) (xj). So the distance between any two different critical points Xj, x 3 - of f\ can 
be estimated from below by 



d(xi, xj) > d(K, e) 



4 K 2 
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Therefore, the number of critical points of f\ does not exceed 



(iii) Suppose that the critical points of fi are x\,--- ,xn,N < N(K,e), and the 
critical values of f\ are ordered increasingly 

fx{xx)<h{x 2 )<...<h{x N ). 

Let g : R — > R be a C fc function satisfying the following conditions 



I 1 ' 


if \t\ 


< 


I o, 


if t 


> 



d{K,e) 

4 ' 
d{K,e) 

d{K,e) ^ |,i ^ d{K,e) 



< ^(t) < 1, if ^ < |t| < , . 
For each i, let Aj : R n — > [0, 1] be defined by Aj(x) = g(||x — Xi\\), and Ci = || || C fc . 
Put r/i = mm(r(K, e), |r^f] ) ■ (The second parameter of min will be used in (v)). 



Let 



JV 

A : R™ R, A(x) = ^CiAi(x), where q = (i - 1) ' /l 



i=i 



Then ||A||c* = maxCj|| \i\\c* < if < lL^£i_ Now consider the approximation of /o : 

/ = fi + A = f + h, where h — I + A. 
We have < ||v|| + ||A|| C * < (e - ^) + rJ ^<e. 

Since Df(x) = iff Dfi(x) = —DX(x), by the definition of A, this equality only 
happens when x G B™ (Jf e) (xj) for some i. But D/i is injective on B™ (K e) (xj) and 

2 2 

D\(x) = DX(x'), when || tjj 0C 2 1 — I 3j *X j 1 1 i SO must be equal to Xj, and then 
Hf(xi) = Hfi(xi) + HX(xi) = Hfo(xi). Thus / is a Morse function having the same 
critical points as fi, and a n (Hf(xi)) > j(K,e) for every critical point Xj. 
Moreover, for any pair of distinct critical points Xj, xj of /, we have 

I/fa) - f(xj)\ = \A(x t ) +*- h(xj) - Cj \ > MK, e) = ^fl. 
We showed that / satisfies (i), (ii) and (iii). 

(iv) Applying the Splitting Lemma 14.21 to / at each of its critical points we get 



5 = 8{K,e) and M{K,e) satisfying (iv). 

(v) First, consider f l = f Q + I. If WDf^x)]] < ±771 = § min(r(if, e), gjffie) )' then 

\\Df (x) — v\\ < r ^' £S> ; and hence, from (i) we have a n (H f\{x)) = a n (Hf (x)) > 
j(K, e). According to the Inverse Mapping Theorem 13.11 Dfx is invertible on a ball 
centered at x with radius p\ = 2 {K+e) 2 ' anc ^ ^ e i ma § e contains the ball centered at 
Dfi(x) with radius p 2 = Pi(-^ + s) > ^1, and hence this ball contains 0. Therefore, 
if ||D/i(a;)|| < \r}\ = | mm(r(K, e), p 2 ), then there exists a critical point Xj of /1 
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such that x G B"(xj). 

Now consider / = / X +A. If \\Df{x)\\ < J^, then ||D/i(x)|| < \\Df{x)\\ + \\DX{x)\\ < 
jrji + \rji = ^rji, and hence x G B™^). Note that Df(xi) = Dfi(xi) = 0. 
Therefore, to get (v) we take 

1 i 2 (K e) -y 2 

n(K,e) = -mm(r(K,e), \) - \ ), and if*(K, e) = mm(5(K, g), 



□ 



Applying the Inverse Mapping Theorem, we get a quantitative version for the 
openness of Morse functions on a ball as follows. 

Theorem 4.5. Let f : B" R be a C k function, k > 2, with ||/||cr* < K. Suppose 
that f is a Morse function with the critical locus £(/) = {xi, ■ ■ ■ , x p } contained in 
B ra and has distinct critical values. Let 

7 = minK( J ff/(x)):xeE(/)} ) 
d = mm{\f(xi) - f(xj)\ : i j and i,j = !,■■■ ,p}, 
p = min( T 2|^2 , ^, d(E(/), <9B n )), and 
= inf{||£>/0i0|| : x G B", d(x, S(/)) > p}. 

Lei e = min(|, |). Then for every C k function f : B™ — >■ R, wni/i ||/ — /Ho* < £ , 
f satisfies the followings: 

(i) If x G B n and ||Z)/(rr)|| < §, then o n (Hf(x)) > |. In particular, / is a Morse 
function. 

(ii) £(/) = {xi, ■ ■ • , x p } C B n , and — Xi\\ < p, for z = 1, • ■ ■ ,p. 

(iii) a n (Hf(xi)) > |, for i = 1, • • • ,p. 

(iv) min{|/(xi) - /(%)| : i ^ j and «, j = 1, • • • ,p} > f . 

Proo/. First note that rj > 0, since £(/) n dP™ = 0. 

Let / : B" -»■ R be a C fe function with ||/- < e. 

Let x G B n , such that < |. Then the definition of e implies 

and hence d(x, £(/)) < p, i.e. ||x — Xi\\ < p. for some i G {1, • • • ,p}. 

Moreover, by ||/||c fc < + £ and the definition of p and e, applying the Mean value 

theorem, we get 

\\Hf(x) - Hf(xt)\\ < \\Hf(x) - Hf( Xl )\\ + \\Hf(xt) - Hf( Xl )\\ 

< (K + e)\\x - x t \\ +e 

< (K + e)p + e 

< (l + p)e + Kp < 2 

Therefore a n (Hf(x)) > ^. This proves (i) and (iii). 

For Xi G £(/), we have < e < |. From (i) and the Inverse Mapping 

Theorem 13 .11 D/ is invertible on a ball centered at Xi with radius 2 (/r+ £ )2 — usk 2 — 
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p, and the image contains a ball centered at Df(xi) with radius 2{k+e) > — 8 > 
\\Df(xi)\\. From these facts and (i), there exists x G £(/) such that \\x — Xi\\ < p. 
Note that, by the definition of p, all critical points of / are contained in B™. 
Moreover, for any two distinct critical points x,y G £(/), \\x — y\\ > 2p. Hence for 
all Xi G £(/) there exists only one Xi G £(/) such that \\x~i — Xi\\ < p, (ii) follows. 
To prove (iv), adding a constant to /, we can assume J(0) = /(0). Then, by the 
Mean value theorem, \f(x) — f(x)\ < e, for all x G B™. So for i — 1, • • • ,p, we have 

\f(xi) - f(xi)\ < \m) - f(Xi)\ + \f(Xi) - f( Xi )\ 
<e + K\\xi - Xi\\ < \ + Kp < \. 

By the triangle inequality, for any two distinct critical points Xi, Xj G £(/), we have 

□ 
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